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THE  MATHEMATICAL  DEVELOPMENT  OF  THE  END-POINT  METHOD 


By  S,  Frankel  and  S.  Goldberg 


ABSTRACT 

The  end-point  method  is  mathematically  developed  and  its  application  to  the  Milne  kernel  studied 
in  detail.  The  general  solution  of  the  Wiener -Hopf  integral  equation  is  first  obtained.  The  Milne 
kernel  appears  in  applying  this  method  to  the  integral  equation  describing  the  diffusion  and  multipli¬ 
cation  of  neutrons  in  multiplying  and  scattering  media.  The  neutrons  are  treated  as  monochromatic, 
isotropically  scattered  and  of  the  same  total  mean  free  path  in  all  materials  involved.  Only  problems 
with  spherical  symmetry  are  treated,  these  being  reducible  to  equivalent  infinite  slab  problems.  So¬ 
lutions  are  obtained  for  tamped  and  untamped  spheres;  in  the  former  case  both  growing  and  decaying 
exponential  asymptotic  solutions  in  the  tamper  are  treated  in  detail.  Appendix  I  treats  the  effects  of 
the  approximations  inherent  in  the  end-point  method  (cf.  LADC  -  79).  Appendix  II  gives  the  solution 
of  the  inhomogeneous  Wiener -Hopf  equation. 


INTRODUCTION 

The  general  development  of  the  end-point  method  and  some  of  its  applications  are  described  in 
LADC  -  79.  It  is  the  purpose  of  this  report  to  supplement  this  general  description  with  an  explicit 
mathematical  development  of  the  end-point  method  and  a  detailed  study  of  its  application  to  the 
Milne  kernel.  This  is  the  kernel  entering  in  the  integral  equation  describing  the  diffusion  and  multi¬ 
plication  of  neutrons  in  multiplying  and  scattering  materials  where  the  neutrons  are  treated  as  mono¬ 
chromatic,  isotropically  scattered,  and  of  the  same  total  mean  free  path  in  all  materials  involved. 

The  end-point  method  of  treatment  of  integral  equations  is  restricted  to  one -dimensional  cases.  This 
essentially  limits  the  method  to  the  treatment  of  problems  in  which  the  materials  involved  and  the 
neutron  distribution  are  both  spherically  symmetric,  these  problems  being  reducible  to  equivalent 
infinite -slab  problems.  In  LADC  -  79  it  was  shown  that  the  end-point  results  may  be  applied  loosely 
to  problems  of  somewhat  more  complicated  geometry  and  give  more  or  less  accurate  approximations 
to  the  truth.  These  applications  depend  primarily  on  loose  analogies  rather  than  mathematical  argu¬ 
ment  and  will  not  be  treated  here. 

Much  of  this  report  will  be,  in  part,  repetition  of  material  treated  in  LADC  79.  Here  the  emphasis 
will  be  primarily  on  the  clear  mathematical  development  of  the  methods  of  application  presented  there. 


Chapter  I 


THE  WIENER-HOPF  METHOD 


The  integral  equation, 


<30 

n(x)  =  J  dx'  n(x*)  K(x  -  x^) 
0 


(1.0) 
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is  known  as  the  equation  of  Wiener  and  Hopf .  With  certain  reasonable  restrictions  on  the  character 
of  K  and  n  this  equation  can  be  solved  exactly.  Before  examining  the  method  of  solving  this  equation 
developed  by  Wiener  and  Hopf,  it  is  useful  to  examine  the  simpler  equation, 


CO 

n(x)  =  j  6x'  n(x')  K(x  -  x') 


Since  this  equation  is  homogeneous,  if  nQ(x)  is  a  solution  then  anQ(x)  also  satisfies  the  equation  for  any 
constant,  a.  Because  of  the  infinite  limits  of  integration  and  the  'Misplacement''  character  of  the 
kernel  (K  depends  only  on  the  difference,  x  -  x'),  nQ(x  -  b)  must  also  be  a  solution.  If  the  solution, 
nQ(x),  is  unique  (except  for  a  multiplicative  factor)  then  nQ(x  -  b)  =  anQ(x)  for  some  a.  Hence  no(x)  = 
ekx^  This  suggests  looking  for  exponential  solutions  of  (1.1), 

oO 

n(x)  =  e*“  =  f  dx'  e^x'  K(x  -  x') 


1  =  6*“  =  /  dx'  6*“'  K(x  -  X') 

~oO 

OO 

=  e^  j  dy  e“ky  K(y) 

-o6 

OO 

j  dy  e-y  K(y)  =  1 


Any  solution  of  this  "characteristic  equation"  gives  a  value  of  k  for  which  e^  satisfies  equation  1.1. 
If  there  is  more  than  one  solution  to  the  characteristic  equation,  then  any  linear  combination  of  the 
exponentials  determined  by  them  will  satisfy  equation  1.1. 

These  considerations  will  be  relevant  to  the  study  of  the  equation  1.0  if  K  decays  rapidly  for 
large  |y|.  If  this  is  the  case,  for  large  x,  equation  1.0  approximates  equation  1.1,  and  it  may  be  ex¬ 
pected  that  with  increasing  x  the  solutions  of  equation  1.0  will  approach  asymptotically  the  exponential 
solutions  of  equation  1.1.  If  this  is  the  case,  the  asymptotic  exponential  part  of  the  solution  of  equa¬ 
tion  1.0  may  be  separated  from  the  remainder  of  the  solution  by  Laplace  or  Fourier  transformation. 
The  use  of  the  Laplace  transform  is  further  suggested  by  the  fact  that  the  left  hand  term  of  equa¬ 
tion  1.2  is  the  Laplace  transform  of  the  kernel. 

Taking  the  Laplace  transform  of  equation  1.1  gives: 

OO  OO  OO 

Jdxe-1“  n(x)  =/  dxe-l“/ dx'n(x')  K(x  -  x') 

-GO  _oO  ~e)0 

OO  OO 

=  f  dx'  n(x')e’^'  f  dye"^K(y) 

-*4o  _oo 

OO  .  OO  . 

j  dx  e"^  n(x)^  j  dy  e"^  K(y)  -  ij  =  0 

This  last  equation  shows  that  the  Laplace  transform  of  n(x)  must  vanish  for  all  values  of  k  which  do 
not  satisfy  the  characteristic  equation  1.2. 

An  application  of  the  same  technique  to  equation  1.0  does  not  lead  immediately  to  a  factored 
equation  because  of  the  finite  lower  limit.  To  get  around  this  difficulty  Wiener  and  Hopf  introduced 
the  following  trick. 

Define  n(x)  =  f(x)  +  g(x) 

where  f(x)  =  0  for  x  <  0 
g(x)  ^  0  for  X  ^  0 

This  permits  writing  equation  1.0  in  the  form 


CO 

f(x)  +  g(x)  =  /  dx'  f(x')  K(x  -  X') 
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Now,  taking  the  Laplace  transform  gives 


Defining 


OO  oo  OO  oo 

J  dx  f(x)  e-^  +  J  dx  g(x)e“‘^  =  J  dx  e"^  Jdx'  f(x')  K{x  -  x') 

-oO  — OO  -«3 

oO  OO 

=  Jdx'  e"^'  f(x')  Jdy  e'^Y  K(y) 

-oo  -OO 

oo 

F(k)  —  Jdx  f(x)  e"^ 

-oo 

oO 

G(K)=  Jdx  g(x)  e"*^ 

—  OO 

K(k)=  JdxK(x)e’^ 


we  have  /  v 

G(k)  -  F(k)  (K(k)-l)  =  F(k)  P(k)  (1.3) 

This  equation  will  hold  for  any  value  of  k  for  which  all  three  integrals  exist.  We  therefore  impose  con¬ 
ditions  on  the  kernel  and  solution  of  equation  1.0,  which  ensure  the  existence  of  a  suitable  region  in  the 
complex  plane  in  which  all  three  integrals  exist.  We  require  that  K(y)  decay  at  least  as  rapidly  as  an 
exponential  for  large  (positive  or  negative)  y. 

K(y)  =  c(e-c|y|),  c  >0.  (1.4) 

Then  K(k)  will  exist  for  -  c<R(k)<  c.  We  further  assume  that 

f(x)  =  c(e^),  d<c  (1.5) 

The  kernels  of  primary  interest  are  symmetric.  For  these,  if  the  “largest''  value  of  c  satisfying 
equation  1.4  is  chosen,  equation  1.5  is  not  a  restrictive  condition,  since  f(x)  must  approach  asymptoti¬ 
cally  an  exponential,  e^,  for  some  k  satisfying  K(k)  =  1  and  therefore  within  the  range  of  convergence 
of  K(k).  The  form  of  equation  1.3  clearly  requires  that  g(x)  decay  (for  large  negative  x)  at  least  as 
fasT as  e^x.  Thus  G(k)  exists  for  all  k  having  R(k)  <  c.  The  three  integrals  will  therefore  all  exist 
throughout  a  vertical  strip  in  the  complex  k-plane  defined  by  d  <  R(k)  <  c. 


Figure  1. 
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Within  this  “common  strip'"  all  three  integrals  are  convergent  and  equation  1.3  must  be  satisfied. 
Outside  this  strip  the  nonconvergent  integrals  will  be  defined  by  analytic  extension  (and  need  not  be 
analytic)  in  such  a  way  that  the  equation  is  still  satisfied. 

Within  and  to  the  right  of  the  common  strip,  F(k)  exists  and  is  analytic,  [it  is  clear  from  its 
definition  that  in  this  range  any  derivative  of  F(k)  exists.]  Similarly  within  and  to  the  left  of  the  strip, 
G(k)  exists  and  is  analytic.  K(k),  hence  also  P(k),  exists  and  is  analytic  within  the  strip  but  may  have 
singularities  on  either  side  of  it.  We  make  the  further  assumption  that  F{k)  and  G(k)  have  no  roots  in 
their  respective  regions  of  analyticity.  (Cf.  Paley  and  Wiener,  Fourier  Transforms,  p.  51).  We 
further  require  that  there  exist  a  sub-strip  within  the  common  strip  within  which  P(k)  has  no  roots. 
[This  must  be  true  if  P(k)  has  only  a  finite  number  of  zeros  in  the  common  strip.  This  will  actually 
be  the  case,  Cf.  Titchmarsh,  Fourier  Integrals,  p.  339.] 

We  have  now  a  sub- strip  within  which  log  P(k)  is  analytic;  within  which,  and  to  the  right,  log  F(k) 
is  analytic;  within  which,  and  to  the  left,  log  G(k)  is  analytic,  and  within  which  the  three  satisfy 

log  P(k)  =  log  G(k)  -  log  F(k) 

This  equation  will  be  satisfied  throughout  the  plane  by  the  analytic  extensions. 

It  is  now  easy  to  find  functions,  F  and  G,  satisfying  this  equation  and  the  analyticity  conditions. 
For  values  of  k  within  the  sub-strip  we  express  log  P(k)  by  means  of  a  Cauchy  integral: 

dk' 

log  P(k)  =  (l/2jri)  log  P(k') 

=  (l/2m)  logP(k') 

+  (l/2«i)  jLk'^k  log 

where  the  contour  of  integration  consists  of  two  vertical  lines  in  the  sub-strip,  one  running  up  to  the 

right  of  k,  the  other  down  to  its  left. 


Figure  2. 
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We  have  now  decomposed  log  P(k)  into  two  parts,  one  certainly  analytic  within  the  strip  and  to  the 
left,  the  other  within  and  to  the  right.  These  may  be  identified  with  log  G(k)  and  -log  F(k),  and  give  a 
solution  to  the  equation  1.0. 


log  F(k)=-^|^^ 


log  P(k)  +  constant 


log  G(k)  =  ^  P{k' )  +  constant 


This  contour  integral  representation  of  log  F(k)  determines  F(k),  hence  also  f(x). 

where  6  is  chosen  to  make  F(k)  regular  along  the  contour..  In  particular,  d  may  be  taken  in  the  sub¬ 
strip.  Since  F(k)  is  analytic  to  the  right  of  the  sub-strip,  the  contour  may  be  translated  to  the  right 
as  far  as  desired.  For  negative  values  of  x  this  may  be  used  to  show  that  f(x)  vanishes. 

If  f(x)  contains  a  term  Aekox  (e.g.,  as  its  asymptotic  solution),  then  its  Laplace  transform,  F(k) 
will  contain  a  corresponding  term. 

oO 

f  dx  e“*^  Ae^ox  = 

Thus  a  pure  exponential  term  in  f(x)  manifests  itself  in  F(k)  as  a  simple  pole,  and  the  coefficients  of 
the  two  may  be  identified.  The  coefficient  of  the  singularity  is  most  easily  determined  by  expanding 
log  F(k)  about  the  singularity. 

log  F(k)  =  -  log(k  -  ko)  +  log  A  +  0(k  -  ko) 

The  asymptotic  solution  will  be  determined  by  all  of  the  singularities  of  F(k)  on  the  imaginary  axis 
and  in  the  right  half -plane.  If  there  are  no  singularities  on  or  to  the  right  of  the  imaginary  axis  the 
solution,  f(x),  will  approach  zero  asymptotically.  A  more  useful  asymptotic  solution  however,  will 
be  that  determined  by  the  first  singularities  to  the  left  of  the  imaginary  axis: 

An  important  special  case  of  this  general  treatment  is  that  for  which  the  kernel,  K(y),  is  sym¬ 
metric  and  for  which  the  characteristic  equation  has  only  a  single  pair  of  conjugate  roots  on  the 
imaginary  axis.  If  these  two  roots  are  at  ±  ik  then  the  solution  will  be  of  the  form 


F(k)  =  Bj^sin  k^  (x  +  x^)  +  h(x)J,  h(x)“^0  as  x- 


Since  the  equation  is  homogeneous,  B  is  undetermined;  Xq,  however,  can  be  evaluated. 


F(k)  =  J  dx  e’^  B  fsin  kQ(x  +  Xq)  +  h(x)1 


2i  W  -  ik. 


e'il'oXo 
k  +  ik_' 


+  2iH(k) 


In  the  neighborhood  of  ±  ik^,  H(k)  is  finite.  We  expand  log  F(k)  near  these  two  poles, 

B 

log  F(ikQ  +  <)  =  log  —  +  i  k^Xg  -  log  €  +  0(e) 


log  F(-  iko  +  e)  =  log  —  -  i  k^Xg  -  log  e  +  0(e) 


limj^log  F(iko+e)  -  log  F(-  ik^^  +e)J  =  log  (-  l)  +  2i  kgX^ 


6 


AECD  -  2056 


log  r{k)  =  log  G(k)  -  log  P(k) 


2»ri  R  k'  -  k 


log  P(k')  -  log  P(k) 


lim  [log  P(iko  +  €)  -  log  P  (-  Iko  +  « )]  =  log  =  log  (*1) 

since  K(y)  is  even,  K(k)  and  P(k)  are  even;  P'(k)  is  odd. 


2  1k(,,„.ijRdk'loeP(k') 


ik  k'  +  iko 


The  two  terms,  log  (-  1),  have  been  neglected  since  the  form  of  the  solution  1.8  is  unchanged  by  the 
addition  of  a  multiple  of  7rto  koXQ.  The  evaluation  of  Xq  completes  the  determination  of  the  asymp¬ 
totic  form  of  the  solution  equation  1.8.  Xq  is  expressed  in  equation  1.9  as  a  single  integral,  which  in 
many  cases  must  be  evaluated  numerically.  To  get  the  complete  solution  requires  two  integrations, 
one  to  evaluate  log  F(k)  by  equation  1.6,  another  to  get  f(x)  by  (1.7). 

Two-Medium  Problems 

A  more  general  problem  that  can  be  treated  by  the  Wiener-Hopf  technique  is 

0  <50 

n(x)  =  J  dx'  K'(x  -  X')  n(x')  +  J  dx'  K  {x  -  x')  n{x'). 

-io  0 

Breaking  up  n(x)  as  before  and  taking  the  Laplace  transform  of  the  resulting  equation  gives 

F(k)  +  G(k)  =  K(k)  F(k)  +  K'(k)  G(k) 
where  the  notation  is  the  same  as  before.  This  may  be  written  as 


G(k)  =  F(k) 


1  ■  K(k)  \  =  F(k)  P(k) 
K'(k)  -  l) 


This  is  now  of  the  same  form  as  equation  1.3.  The  rest  of  the  treatment  proceeds  in  the  same  way. 
With  this  more  complicated  form  for  P(k)  there  may  be  a  greater  number  of  singularities  of  log  P(k), 
leading  to  a  larger  number  of  independent  solutions.  In  particular  it  is  no  longer  necessary  to  re¬ 
quire  that  g(x)  decay  exponentially  away  from  the  boundary. 

An  important  special  case  of  this  two-medium  problem  is  that  for  which  K(y)  and  K'(y)  differ  only 
by  a  multiplicative  factor.  This  case  will  be  treated  extensively  in  the  second  chapter. 

The  Wiener-Hopf  technique  may  be  further  extended  to  permit  the  solution  of  inhomogeneous 
displacement  integral  equations.  This  method  is  outlined  in  Appendix  II. 


Chapter  11 


APPLICATION  TO  NEUTRON  PROBLEMS 


In  this  chapter  we  treat  the  applications  of  the  Wiener-Hopf  method  (combined  with  some  approxi¬ 
mations)  to  problems  concerning  the  spatial  distribution  and  time  dependence  of  neutrons  in  spheres 
of  multiplying  and  scattering  materials.  It  will  be  shown  that  such  problems,  with  suitable  physical 
approximations,  can  be  represented  by  integral  equations  closely  analogous  to  the  Wiener-Hopf  equa¬ 
tion.  By  making  suitable  mathematical  approximations  (the  ''end-point  method^ 0  fairly  accurate 
solutions  to  these  equations  can  be  gotten  from  the  corresponding  Wiener-Hopf  solutions. 
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We  make  the  following  physical  approximations: 

A)  We  consider  only  one  neutron  velocity;  hence  for  each  material  only  one  value  for  each  cross 
section. 

B)  We  treat  all  collision  processes  as  isotropic.  (Anisotropy  of  elastic  scattering  can  be  treated 
to  a  limited  extent.  It  can  be  shown  that  if  this  anisotropy  is  neglected  and  the  transport  average  used 
for  the  elastic  scattering  cross  section  quite  accurate  results  will  be  obtained.  Cf.  LADC  -  79  and 

MT  -  26.) 

C)  The  total  mean  free  path  will  be  taken  to  be  the  same  for  all  materials  involved. 

D)  The  neutron  distribution  will  be  treated  as  a  continuum.  It  will  be  taken  to  be  spherically 
symmetric  and  of  stable  spatial  distribution.  These  three  conditions  will  certainly  be  good  approxi¬ 
mations  if  the  neutron  distribution  has  lived  through  many  generations  and  consists  of  a  sufficient 
number  of  neutrons  to  make  statistical  fluctuation  negligible. 

We  adopt  the  following  notation: 

Of  is  the  fission  probability  per  unit  path  length.  (It  is  therefore  the  product  of  the  fission  cross 
section  and  the  number  of  nuclei  per  unit  volume.)  Similarly, 

(Tg  is  the  scattering  probability  per  unit  path  length. 

(Tg,  is  the  absorption  probability  per  unit  path  length. 

(T  =  +  C7g  +  (7g 

V  is  the  mean  number  of  neutrons  emerging  from  a  fission  process. 

VOf  +  <To 

F  =  1  +  f  = - - - is  therefore  the  mean  number  ux  neutrons  emerging  from  a  collision. 

V  is  the  neutron  velocity. 

n(r,  t)  is  the  neutron  density  at  point  r  at  time  t. 

We  express  the  neutron  density  at  {r,t)  as  an  integral  over  all  points  at  which  these  neutrons  may 
have  suffered  their  last  collisions. 

vn(r,t)=/dr'  <tv  F(r' )  n(r',  t  -  r'  )2  ^  '  (2.1) 

We  look  for  solutions  of  the  form 

n(r ,  t)  =  n(r )  e^ 

The  integral  equation  2.1,  then  takes  the  form: 

We  now  rescale  r,  taking  as  the  unit  of  length  the  mean  attenuation  distance,  l/(<r  +7  q/v). 

x  =  t(c  +  Vq/v) 


“(2)  =  1  ^r\/cv 


1  +ro/ffV  j  -  "  '  4n[x  -  X')2 

Defining  7-7  o/tfv  gives  the  three-dimensional  integral  equation 

1  r . 


n(x)  =  Jdx'  F(x')n(x') 


4n(x  -  x')2 


'  If  we  now  introduce  polar  coordinates,  x'  =  (r',  0')? 


taking  the  point  x  on  the  polar  axis  we  may  make  use  of  the  assumed  spherical  symmetry  of  n(x')  to 
reduce  equation  2.2  to  an  equation  in  one  dimension. 

-  -  (r2  +  r'2  -  2rr'  cos  0)1/2 

n(r)  =  — —  r^  2  dj.  ^  F{rO  n(r ')  I  d0sin0d0"~7^ - — r — - — — — 

1 +y  -  '  '  '  '  Jj  ^  47r(r'^  +  r'2  -  2rr'  cos^ ) 
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^  [E(|r-r'|)-E(r  +  r')] 


rn(r)  =  -r77^ - ;  f  dr'  F(r' )  r'  n(r' )  [E(|r  -  r'l )  -  E(r  +  r' )]  (2,3) 

2(1  +  y ) 

0  , 

If  we  now  define  u(r)  =  r  n(r)  and  treat  u(r)  as  an  odd  function,  and  F(r)  as  an  even  function  of  r  [no 

meaning  has  previously  been  assigned  to  negative  values  of  r  or  to  the  corresponding  n(r)  and  F(r)]  , 
we  may  write  equation  2.3  in  the  form: 


u(r)  =  2(1^+ 7)  F(r')  u(r')  E(lr  -  r'l)  (2.4) 

If  instead  of  assuming  the  material  and  neutron  distribution  spherically  symmetric,  we  take  both  as 
functions  of  only  one  Cartesian  coordinate,  z,  equation  2.2  may  be  reduced  to  an  equation  in  one 
dimension  as  follows: 

^  ^  [(z  -  z')2  +  (y  -  y')2  +  (x  -  x')2]^ 

n(z)=j^/dz'  F(z')n(z')  ffdx'dy'  .  ,.)2  +  (y  -  y')2  .  (x  -  x')2J . 

27r  «  g.i 

"TTr  /dz'  F(z')  n(z')  J  d<A  Jpdp^2 

o  0 

where  1^  =  (z  -  z' ^  ^ 

n(z)  "  2(lVr~)  F(z')  n(z')  E(|z  -  z'l)  (2.5) 

A  comparison  of  equations  2,4  and  2,5  shows  that  the  sphere  problem  2.4  may  be  identified  with  a  slab 
problem  2.5  in  which  the  distribution  of  materials  F(z)  across  the  slab  is  the  same  as  that  along  a 
diameter  of  the  sphere.  Any  odd  solution  of  the  slab  problem,  n(z),  may  be  identified  with  the  quantity 
u{r)  in  the  sphere  problem  and  conversely.  The  “fundamental  mode'’  of  the  sphere  for  which  n(r)  is 
everywhere  positive  corresponds  to  the  “first  harmonic"  of  the  slab  in  which  the  neutron  density  takes 
on  apparently  meaningless  negative  values.  For  this  reason,  and  because  higher  modes  may  be  super¬ 
imposed  on  the  fundamental,  we  will  treat  the  neutron  density,  n(z),  as  a  real  quantity  which  may  have 
either  sign. 

For  a  tamped  sphere  of  core  radius  a,  outer  tamper  radius  b,  and  mean  attenuation  distances,  the 
integral  equation  2.4  takes  the  form 

-a 

u(r)  =  I  ^  f  dr'  u(r')-|-E  (Ir  -  r'l) 
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+  j  dr'  u(r')iE(|r-r'|) 

+  J  u(r')  i-E  (|r  -  r'l) 

a 

where  and  ft  are  the  values  of  f  in  core  and  tamper  respectively.  This  equation  differs  from  the 
Wiener~Hopf  equation  in  having  four  boundaries  instead  of  one  (or  two  for  an  untamped  sphere).  With 
more  than  one  boundary  no  exact  solution  is  known.  We  therefore  resort  to  an  approximation,  namely 
to  treat  the  behaviour  of  the  solution  near  each  boundary  as  if  no  other  boundaries  existed.  It  was 
shown  in  the  first  chapter  that  the  solution  of  the  one-bovindary  problem  approaches,  at  large  distances 
from  the  boundary,  a  solution  of  the  problem  with  infinite  limits.  It  is  reasonable  to  expect  that  the 
solution  of  a  two-boundary  problem  in  which  the  boundaries  are  very  far  apart  will  behave  in  some 
middle  region  as  a  solution  of  the  infinite -limits  equation.  If  this  is  the  case,  we  have  only  to  combine 
two  one-boundary  solutions  in  such  a  way  that  their  asymptotic  components  coincide.  In  a  many¬ 
boundary  problem,  e.g.,  the  tamped  sphere,  we  apply  this  recipe  in  each  region.  This  approximation 
method,  the  “end-point  method'^,  would  seem,  from  the  above  argument,  reasonably  accurate  only  if 
the  distances  between  boundaries  are  many  mean  attenuation  distances.  It  is  shown  in  Appendix  I 
that  the  limit  of  reasonable  accuracy  is  actually  a  few  tenths  of  a  mean  attenuation  distance.  There 
is  therefore  good  reason  to  believe  that  for  sizes  larger  than  that,  the  end-point  method  is  sufficiently 
accurate. 

In  order  to  apply  the  end-point  method  we  must  first  study  the  one -boundary  problem  with  the 
“Milne  kernel’ ^ 

K(y)  =  c-|.E(|y|) 

This  kernel  with  c  =  1  occurs  in  ^^he  equation  of  E.  A.  Milne”  describing  the  flow  of  radiation  through 
the  outermost  layers  of  a  star.  We  will,  however,  refer  to  it  as  the  “Milne  kernel”  for  all  positive 
values  of  c.  The  general  equation  we  have  to  study  is  then 

n(x)  =  c'  J  dx'  n{xO  "^Edx  -  X' I)  +  c  j*  dx'  n(x' )  ^  E  (|x  -  x'|) 
c  =  (l-ff)/{l+'y). 

Several  cases  arise.  For  a  free  surface,  either  the  outer  surface  of  a  tamper  or  the  surface  of  an  un¬ 
tamped  sphere,  we  take  c'  =  0.  For  an  interface  we  take  both  c  and  c'  positive.  For  the  core  material, 
c  must  be  greater  than  1  (f>  V);  in  the  tamper,  c  -  1  may  be  of  either  sign. 

We  first  treat  the  free -surface  case. 

OO 

n(x)  =  c  J  dx' n(x')-|-E  (|x  -  x'i ) 

The  characteristic  equation  is 


1 
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If  c  <  1  we  have  two  real  roots,  ±  such  that  c  =  k/tanh  If  c>  1  we  have  two  imaginary  roots,  ± 
i  ko,  such  that  c  =  kQ/tan“*lko,  In  either  case  it  can  be  shown  that  the  characteristic  equation  has  only 
two  roots.  In  the  latter  case  the  asymptotic  solution  is  a  sinusoidal  function  of  k^x,  in  the  former,  a 
hyperbolic  fimction.  We  will  represent  the  phase  of  the  asymptotic  solution  by  the  ' ‘extrapolated  end¬ 
point, Xq,  such  that  the  asymptotic  solution  is  the  sine  or  hyperbolic  sine  of  kQ(x  +  Xq).  We  now 
follow  through,  explicitly,  the  method  of  solution  outlined  in  Chapter  1. 


oo  ^ 

n(x)  =  f(x)  +  g(x)  =  c  I  dx' f(x')'2  E  (|x  -  x'l) 


f(x)  =  o  for  x  <  o 
g(x)  =  0  for  X  ^  o 


CO  Ir  W 

F(k)  +  G(k)  =  J  dx  n(x)  e  =  J  dx  e  J  dx '  f(x' )  |  E(|x  -  x' | ) 


j  dx'f(x0e^  J  dye~^|E(|y|) 


=  F(k)  ^  log 


G(k)  =  F{k)|—  log  (yrij  “  ^ j- 

P(k)  has  singularities  only  at  ±  1.  These  singularities  are  branch  points  so  that  to  make  the  fimction 
explicit  we  introduce  cuts  lying  along  the  real  axis  from  to  -1  and  from  +1  to  +«?.  We  treat  first  the 
case  c>  1.  The  two  roots  of  P(k)  are  then  pure  imaginary,  ±  ik^.  The  singularities  of  log  P(k)  are±l 
andlikp.  We  look  for  a  log  F(k),  analytic  to  the  right  of  the  imaginary  axis  [corresponding  to  the 
sinusoidal  asymptotic  solution,  f(x)],  and  a  log  G(k),  analytic  to  the  left  of  +1  [corresponding  to  a  g(x) 
decaying  somewhat  faster  than  e-xj,  and  satisfying 

log  P(k)  =  log  G(k)  -  log  F(k)  (2.6) 

The  ‘‘sub-strip''  in  which  all  three  of  these  quantities  are  analytic  is  0<R(k)  <  1.  We  therefore  break 
up  log  P(k)  by  means  of  a  Cauchy  integral  along  a  contour  running  up  and  down  in  this  strip  and  en¬ 
closing  k,  and  (except  for  a  common  constant)  identify  log  G(k)  and  -log  F(k)  with  the  two  parts  of  the 
integral. 


log  PR(k) 


R  k'  -  k 


-  log  P(k')  =  log  G(k)  +  constant, 


log  PL(k) 


f  dk '  log  P(k' 
L  k'  -  k 


log  F(k)  +  constant. 


We  simplify  log  PR(k)  by  deforming  the  right  contour  to  enclose  the  right-hand  cut. 
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Here  the  tan“l  rises  from  0  at  k'  =  1  to  tt  at  k'  =  +  «  (as  indicated  by  the  bracketed  expressions). 
Substituting  k '  =  1/s, 

Where  _  i  /  »r/2  \  Tc  =  ir  s  =  0 

^  “  \tanh-ls  -  1/cs/  =  0  s  =  1 

log  PH(k)  -  I  J‘  f  T,  ♦  S  To  (2.7) 


Here  and  throughout  this  treatment  we  encounter  logarithmically  infinite  constants.  A  slight  modifi¬ 
cation  of  our  procedure  [to  make  P(k)  *-^1  as  I  k|— ooj  suffices  to  avoid  this  embarrassment.  The 
present  treatment  is  somewhat  simpler,  though  formally  less  rigorous. 
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We  simplify  log  PL(k)  by  a  corresponding  deformation  of  the  left  contour. 


Letting  s  =  ”  we  have 
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-  log  PL(k)  =  -  1  /  [  2»r  +  tan-1  j  j^tan’l  =  -  Tc  =  -»r— o] 

7T 


k2  +  ko2 
(1  +  k)2 

/  ds  1  r 


■■*/  f'ij  ShTS 

0  0 

log  PL(k)  =  2  I  y  i  J  f  Tc  -  log  (k2  +  ko2)  +  I  I  ^ 


Combining  those  two  expressions,  2.7  and  2.8,  with 


log  P  (k)  =  log  log  -  l)  =  log  PR(k)  -  log  Pl  (k)  (2.6) 

+  log  (k2  .  ko2)  +  ^  ^  1  -  k2s5  '^c(2.9) 


Taking  the  limit  as  k--i-0  we  get 


1  r  ds  V  1 , 

J  -(Tc-»r)=^log 


and  equation  2.9  becomes 


^  4  sds  _  ^  ds  k  4  ds 

*■  1  -  k2s2  ”n  l-ks  ‘^"wj  1  +  ks  ' 

4  o  o  o 


Dividing  by  k^  and  again  letting  k— ►O, 


1 ,  /  k2  +  ko2  ^  1  , 


2i2"6(^^) 


c  ,  l±i  , 

2k  1  -  k  ■  ^ 


r  1  ds  1  rig 

We  now  subtract  the  (infinite)  constant,  2  ” —  3  ~  ~  from  log  PR(k)  and  log  PL(k)  to 

JOS  ^  J  o  ® 

give  log  G(k)  and  log  F{k). 


log  F(k)  =  -  log(k^  +  kg)  +  I  f  Tp  +  log  B; 
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logG(k)  =  |  [  ^  (Tc-»r)H-|  f  Tc  +  logB, 


k  r  ds 


k  r  ds  ^  - 

”  tr  J  1  -  ks  ^  ko^ 


B(c  ^  1) 


We  now  determine  Xo  and  the  value  of  B  required  to  give  the  asymptotic  sine  wave  in  f(x)  unit  ampli¬ 
tude. 


f(x)  =  sin  ko(x  +  Xq)  +  h(x) 


h(x)— *-0  as  X— *-+ 


F(h.)  -  +  H(k)  =  ksin  kpxp  +  ko  cos  kpXQ  ,  , 

log  F(iko  +  « )  =  -  log(2i)  +  ikoXo  -  log  e  +  0(<) 

log  F(-iko  +  «  )  =  -  log(-  21)  -  ikoXo  -  log  €  +  0(e) 

lim  [log  F(lko  +  e )  -  log  F(-  Iko  +  e )]  =  log  (-  1)  +  21koXo 
€-•0 


1 

=  ^  J  ds  Tc  (z — ^  +  ■; — +  log  (-  1) 

7t  ^0  ^  \1  +  ikos  1  -  ikQsy 

1 

1  f  ds  „ 

^  J  TTkJ^  ^ 


ds  „ 

+  k  2s2 


Now  adding  the  two  values  of  log  F  gives 


log  F(iko  +  €  )  +  log  F(iko  +  < )  =  -  2  log  (2«)  +  O(^), 

=  2  log  (2koel  ds  Te(Y^  •  +  2  log  B  +  0(e) 

2krt2  ?  Cl  rfc 

=  -  2  log  (2ko<)  +  -^  J  Tc  +  2  log  B  +  0(e). 

.  _  .  .  ko2  ,1  s  ds  „ 

log  B  =  log  ko  =  —  “k^sZ  Tc 


This  integral  may  be  evaluated  by  allowing  k  to  approach  ikp  in  equation  2.11: 
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k„2  f  s  ds  _  ,,  r  1  ,  /'2iko«\  1  .  -  c  -  1 


1  .  2(c  -  1) 

2^°^  c' 

■  1  +  ko^ 

®  ®  2  2(c  -  1) 


logF(k)  =  ^  j  Tc  -  log  (k2  ko2)  ^ylog  ^ 


p/w - i5a - \fTZs/K±^^  ff/ol-i-ks 

k2  +  ko2  V  2(c  -  1)  e 


k  1  ds 

1  -  c/(l  +  k„2)  IT  j  2  +  ks 
— rr^ — rr“^e  Jn 


k  sin  ko  Xq  *  ko  cos  kQ  Xq 


Wfc  can  evaluate  H{o),  the  total  area  of  h{x),  and  *'mean  length*', 


H(o)  = 


I  1/ 1  -  c/(l  +  k 


~  U^nMr.2  I  *'0*0  ■  *^1 


H(o)  H(o)ko 
Making  use  of  the  formula 


yi  -f/f.ii.ij'o!) 

0  V  2(C  -  1) 


-if 

n  Jo 


n(o)  llm  k  r  dx  n(x)  e  11m  k  F(k). 

k-«  i  k-^* 


we  get 


k— 


»>»'  -  »*  ^  ■  1^331^. 

We  can  derive  an  expression  for  h(x)  suitable  for  numerical  evaluation  as  follows: 


1  i®®  +  0 

h(x)  =  ^  f  dk  e™  H(k),  0  <  8  <1 

2  1  J  4_  ,  « 
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H{k)  is  not  singular  at  +  iko.  The  bracketed  expression  vanishes,  thus  the  contour  may  be  deformed 
to  lie  along  the  left  cut.  Only  the  integral 


^  I  T  -  ~  r  2^  r  s  ds 

^  Jo  1  +  ks  ‘  ^  1  -  ks  ‘  ^  Jo  1  .  k2s2 


-  f- 

^  J  1 


Tc  “  log 


c  ,  1  +  k  . 

ko2  2k  1  -  k  ~  ^ 
k2  +  ko2  c  -  1 


is  double-valued  across  the  cut.  Thus  only  the  first  term  in  H(k)  contributes. 


kx  ko  i/ 1  -  c/(l  +  ko2)  tt  J  i  -  ks  (c  -  l)(k2  +  kn^) 

k2  +  ko2  V  2(c  -  1)  ^  ko^ 


{h(x)  is  negative  for  all  x). 

If  c  <  1  the  roots  of  the  characteristic  equation  are  ±  kj,  where  c  =  kj^/tanh'^kj^.  The  contours  must  now 
be  taken  as  shown  in  Figure  5. 


Figure  5. 
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Proceeding  in  the  same  way  as  for  c>  1  we  get  the  analogous  results: 

n(x)  =  sinh  k^{x  +  Xq)  +  h(x) 


1  f  V  i-  ,  A  - 

^  J  -  (Tc-,r)  =  2log^ 


1  ,  1“C 


f  sds  _  1  ,  kp^  -  .  1  -  c 

TT  J  1  -  k2s2  c  “  ~2  ki^  c  .  1  +  k 


1  ,  ko2  -  k2  .  1  -  c 


r _ _ T 

0  ^  Jo  1  “  kj2s2  < 


tanh"l  s  -  1/cs 


r^r[/rs’  [tan-i=^o] 

k  /'I  ds 


p,l,\  -  —k] _ y'c/d  -  ki^)  ^  F  / 

F(k)-,,2-ki2V  2(1 -c)  ' 


1  +  ks 


“W'-Fi 


1 

H(0)ki' 


c/(l  -ki2)  -  1 

2(1  -  c) 


_  k  fl  ds 

-  r°^  ^1,  '  Jo  1  +  ks  ^ 

j  c  \  j  k  dk  e  ^ 


Combining ;these  hyperbolic  results  (c  <  1)  with  the  elliptic  results  (c  >  1)  previously  obtained  shows  the 
character  of  the  solution  and  its  numerically  identifiable  features  to  be  continuous  (as  a  function  of  c) 
across  the  parabolic  (c  =  1)  boundary  case. 

We  now  treat  the  two-medium  case,  distinguishing  the  two  materials  (e.g.,  active  material  and 
tamper)  only  by  their  different  values  of  c.  Here  four  cases  arise  as  the  two  c  values  are  less  than  or 
greater  than  1.  We  treat  explicitly  only  the  case:  c  >  1,  c'  <  1.  The  extension  to  other  cases  will  then 
be  obvious.  Because  of  the  applicability  of  the  solution  to  the  simple  tamped  sphere  we  refer  to  the 
one  region,  c>  1,  x  >0,  as  ^^the  core^',  and  to  the  other,  c  <  1,  x  <  0,  as  “the  tamper'\  We  find  two 
pertinent  solutions,  one  belonging  to  a  growing  and  the  other  to  a  decaying  exponential  asymptotic  so¬ 
lution  in  the  tamper.  For  the  problem  of  the  infinitely  tamped  sphere  only  the  decaying  solution  will 
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figure  (decaying  as  one  moves  away  from  the  interface  into  the  tamper).  However,  the  ‘^asymptotic 
solution”  for  a  finite  tamper  will  be  a  linear  combination  of  the  two  solutions.  The  integral  equation  is: 

n(x)  =  c'  ^  n{x')-|E(jx  -  x'| )  +  c  J  n(x')  -ifiOx  -  x'l)  (2.14) 


We  use  the  same  notation  as  before: 


n(x)  =  f(x)  +  g(x) 
f(x)  =  0,  X  <  0 
g(x)  =  0,  x^O 

F(k)  =  f  dx  f(x)e  ^ 

*'-0O 

G(k)=  f  dx  g(x)e 

.Leo 

K(k)=  J  dxi-E(ixl)  ^  log|-^ 

F(k)  +  G(k)  =  I  dx  n(x)e*^ 

•'-•O 

-00  -OO 

=  J  dx  e"*™  j  ^  dx'  -jE  (|x  -  x'l  )  [^c'  g(x')  +  c  f(x')J 
=  J*  dy  e-ky|E(|y|)  dx’e'*“  j^c'g(x')  +  c  f(x‘)J 

=  ^  log  [o'  G(k)  +  c  F(k)] 

c  ,  1  +  k 

^  log  yzy  ■  ^ 

G(k)  =  F(k) -  ^  f  =  F(k)  P(k) 

The  singularities  of  log  P(k)  now  lie  at: 

±  1  (branch  points) 

±  iko  [roots  of  P(k), 

-  kf  [poles  of  P(k)^^^ri^  =c'] 

F(k)[and  we  assume  also  log  F(k)  J  must  be  analytic  for  R(k)  >  0 
G(k)^and  we  assume  also  log  G(k)  J  must  be  analytic  for 

R(K)  <  +  k]^  for  “decaying  solution”,  i.e.,  g(x)  =  O(e^l^) 
or  R(k)  <  -  k^  for  “growing  solution”,  i.e.,  g(x)  =  0(e''^l^) 
log  P(k)  is  analytic  for  -  1  <  R(k)  <  +  1,  except  at  ±  iko,  i 
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We  treat  first  the  decaying  solution.  As  before  we  identify  log  F(k)  and  log  G(k)  with  the  left  and  right 
integrals  (again  excepting  a  constant). 

log  PR(k)  =  ^  J  log  P(k' )  =  log  G(k)  +  const, 

log  PL(k)  =  "  ^  J  P(k')  =  log  F(k)  +  const. 

L 

We  deform  the  contours  as  follows: 


Figure  8. 


making  use  of  the  previous  evaluation  of  the  first  term. 
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Figure  9. 


The  last  integral  is  now  equivalent  to  that  evaluated  in  equation  2.15  (and  is  identical  with  the  right- 
contour  integral  occurring  in  the  one -medium  problem  for  c  <  1). 


ds 

Jjj  s(l  -  ks) 


,  1  ^ 
i  r  T 

■n  Jo  s(l  -  ks)  ^ 


k  /  ds 


‘'O 


^0  o 


We  choose  the  constant  to  make 


1  r  ds  /  r  '  ^  ds 

log  G(k)  =  log  PR(k)  +  log  B  -  —  (T^.  -  0  -  I  ” 

‘'o  o 

Evaluating  the  left-contour  integral  gives 

-  log  PL(k)  =  ^  /  Hk^  ■  1  -k'  ) 


1  r  ds 


=  -2  f  +log(k2  +  ko2)+-  r  ^ 


^ks) 


AECD  -  2056 


21 


Figure  11, 


1  r  dk  .  /  c'  1  +  k'  ,\  If  dk"  f  c  .  1  k"  \ 

Ikl  k'  -  k  (2k'  1  -  k'  ’  ”  2«i  +  k  2k"  1  -  k  "  *  7  ’ 


=  lf  ^ 

If  J  s(l  + 


s(l  +  ks) 


.logPL(k)  =  -2  J  ^+log(k2  +  ko2)+i  J  ^  r^)  -i  /i(Ste)Tc 

0  00  o 


1  l/k, 

ft  1 


o  o 


log  F(k)  =  log  PL(k)  +  log  B  -  i  J 


(Tp-Tc) 


k  r  ds  _  ,  ,  (kl  +  k)B  2  f 

"irJl  +  ks  *^  c')'*'®?  ki(k2  +  ko^)  It  J 


(Tc  -  Tc' )  +  2 


/  ? 


I  {?  [("-Tc')-(’f-Tc)]=-logjl:^ 


=  I  J  irli  -  Tc^)  +  log  -  108  •  !?t)  "  •'I' 


_k  f  ds  «  X  ,  Wm±Mi_LO 

-  TT  J  1  +  ks  1^0  -  Tc')  +  log  +  ko2){c  -  1) 


Bko2(ki  +  k)(l  -C-) 
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We  again  determine  xq  and  the  value  of  B  required  to  make  the  asymptotic  sine  solution  of  imit  ampli¬ 
tude. 

f(x)  =  sin  ko(x  +  Xq)  +  h(x),  x>  0,  h(x)-*-0  as  x-*-+«‘>  (2.17) 

Ik-iko  '  k+iko 


lira  [log  F(iko  +  «)  -  log  F-iko  +  «)]  =  log  (-1)  +  2ikoXo 


=  ^  I  1  (’’c  -Tc')  +  log(^|^)  +  log 


ki  +  iko 
ki  -  iko 


If  ds  X  1  .  -1  ko  1  ,  -1  kn 

=0  =  ;^  Jm^(Tc-Tc.)+i^tan  ^  =  xj  +  -  tan  1 


lim  [log  F(iko  +  « )  +  log  F(-iko  +  «)  -  2  log  <]  =  -  2  log  2 


2ko2f  _j 

TT  J„  1  + 


s  ds  X  «  .  Bko2(l  -  c')  ,  ki2  + 1 

0 17^  (Tc  -  Tc<)  .  2  log  --g^Tir-  - 


The  first  term  may  be  evaluated  by  the  use  of  equation  2.11  and  equation  2.13. 


’‘^o^  r  sds  ,  2ikoe(c  -  1) 

»  J  1  +  ko2s2  '  c  '  ^  ^21.  1. _ £. 


2(c  >  l)ki2(l  -  c  Vc) _ 


.  «  .  ki  (c  -  2)  1  ,  ki2  +  kpg  1  ,  2(c  -  1)  ki2(i  -  c'/c) 

ko2(l  -  c') ’  2  ko2  2  [l  -  c/(l  +  ko2)]  (ki2  +  ko2)(l  -  cO 


1.  (c  -  dFi -c/(l  +  kn2)' 

2  2ko2(l  -  c*)(l  -  c'/c) 


log  F(h)  = 


_k  f  ds 

=  w  J  1  +  i 


(Tc  -  Tc-)+ »  log 


1  ko^d  -  c')[l  -  c/(l  +  ko2). 


2  ®  2ki2(c  -  1)(1  -  c'/c) 


/_k+J^N 

\k2  +  ]s^2j 
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F(k)  =  ^  k  +  ki  JiT-  C-)  [l  -  c/(l  +  k^2)]  n  f  1  -  Tc') 

ki  k2  +  ko2  V  2(c-l)(l-c'/c)  ®  { 

H(k)  =  F(k)  “  ^  sin  kpxp  +  ko  cos  kpXQ 
k2  +  ko^ 


I - r - T  r  ds 

^  i..vl/U’C0[l~c/(l-Kkn2)]  J  1  +  ks  -  T, 


k2  +  ko2  2(^-T)(1-cVc) 


-  k  sin  koXo  “  kQ  COS  k^Xo 


H(0)  =  ^  Ll-:  c')ll  -c/(l  +  kn2)] 

'  '  ko  2(c  -  1)(1  -  c'/c)  ^*0 


=  /<^(Tc-Te0]  -^sinkoxo 


iOi) _ L_  ./(l  -cQ  I  1  -  c/(l  +  kn2)  I  /  J_  ^  1 

H{0)  "  H(0)  ko  V  2(c  -  1)(1  -  c'/c)  Ui  n 


U1  Jds(T^-To.)j  " 


sinkoXo 


n(0)  =  Um  kF(k)  =  lim  k  (k  +  k  -  c')[l  -  c/(l  +  kn2)]  w  V) 


k^«  k2  +  k^2  ki  ^  1'  V  2(c  -  1)(1  -  C  '/c) 


»  I  1  As^'^c  ■  "^c')  iJ  T  -")  -sl^  ~(Tc'  -w) 


(c  -  1)  ki2 
^2(1  -  c') 


(using  (2.10)  and  (2.12)). 


1  -  c/(l  +  ko2) 
2(1  -  c'/c) 


h(x)  =  ^  f  H(k) 


^  i 

ITTI  J  ri> 


h)n”  /  1  +  ks^'^C-Tc') 


j  dke^(k  + 

Ijj’  k2  +  kQ2 


Ce  •'o 
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^ 

-1 

1 

Figure  12. 


where  C  = 


^  l/Ci  -  c/(l Tko^)]  (1  -  c') 
ki  V  2(c  -  1)(1  -  c  7c) 


k(x)  =  'f  dk  e“  C  ^  "S  e  •'o  ‘  j  - 

2«i  ko2(kj  -  k)(l  -  c')  I  £ 


(Tc  -  Tc') 


dk  e*“  C 


- 


S2_  „  r  1  <1=^  *  ~  2k  1  -  k  ) 


k  r  ds  S  .  C'  A  1  4-  k 

s -^v 

°  c  /,  1  +  k  . )  , 


Replacing  k  by  -k  gives 

i-w-sr  /*'■ 


1  - 

-( 

2k' 

1  4-  k 

1  -  k 

4-  Tfi^ 

5' 

(log 

1  4-  k 

1  -  k 

4-Tri 

-kx  ki  ;  /[I  -  c/(l  +  ko2)]  (e  -  1)  ^  "  tt  j  1  +  ks  ) 


ko(ki  +  k)V  2(1  -  c')  (1  -  cVc) 


k  ?  ds 


where 


/.  c'  k  +  1\  ,  /  c  1^  k  4-  1  \ 

d  c  V  '  2k  k  -  ih  (  2k  k  -  1  ~  7  _  *[i _ c  -  c' _ 

c  /c  .  k+1  ,\2  /  cjr\2  "  '  k  /  c  k  +  1  7  2  ,  c»r\ 


>F - 

, ,  V  kA/Ll  -  c/(l  +  kn2)]  (c  -  1)(1  -  cVc)  f  kdk  _ e _ ° _ ^ 

2(1 -c')  •{  k  +  ki  j ^  2  ^ 


Now  returning  to  G(k)  .  1  . 


k  i  ds  .  rr.  \  1  ki  1  (c  -  1)  [l  "C/(l4-kf^2)~| 

7T  /  1  -  ks  “  Tc' )  +  log  .  1.  +  2  2ko^(l  -  c0(l-c'/c) 
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A  check  of  this  expression  is  afforded  by  evaluating 

g(-  c)  =  lim  -  k  G(k)  =  n(0), 

k— “oo 


(ef.  equation  2.20), 


G{k)  =  r  dx  e  ^  g(x)  =  f  dx  e  [ab  +  j(x)]  , 

•f- to  •l-«o 

where  j{x)  =  0(e  as  x—- 
G(k)  =  ^  +  Jlk),  J(k^)  is  finite. 

log  G(ki  +  €)  =  log  (  7)  +  0(c) 

(rkl\3  f  ds  (c^l)[l-c/(lH-kf,2)] 

\  s  /  TT  J  1  -  k^s  ^c  ^  +  2  2ko^(l-c0(l-c  Vc) 


r - - -  M  r  ds  ,  . 

A  \  +  ko2)]  TT  J  1  -  k^s  ^  ^ 

^  ko  \  2(1  -  c0(l  -  c'/c) 


1  1  o 

ki  r  ds  „  X  ki  r  ds  ,  ki^  s  ds  , 

I  1  -  kis  "TJq  1  -  ki2s2  (Tc  -  Tc')  +-^  1  .  j^,2s2  -  TcO. 


TT  1  ~  ki2s2 


The  first  term  will  be  called  kix2  by  analogy  with  the  xj  introduced  in  equation  2.18,  the  second  can 
be  evaluated  by  the  use  of  equation  2.11  and  2,13. 


so  that 


^  It^ATc-Ic) 


2kn2(c/c>  -  1)(1  -  C') 


(ki2  +  ko^)(c  -  1)  [c  /  (1  -  ki2)  -  i] 


[1  -  c/(l  +  kn2)]  kix2 


ki^  +  k^^  c'l^c  V(1  -  kj^2) -1^ 

ki  Vc[l  -  cAl  +  kn2^ 

■Vki2  +  ko2  yc'  ~c'/(l  -  kj^)  -  1' 


ki{x  +  xi) 

,e  ^  +](x) 


J(k)  =  G(k)  - 


kx  -  k 


ki  ./7-  1)  [1  -  C/(1  +  kTT 
ko(kx-k)V  2  (1  -  c')(l  -  c'/c) 


26 


AECD  -  2056 


icn 

'W. 

-ioc 


r' 


ko2(ki2  -  k2)(l  -  c')(|^log  l) 

(k2  +  ko2)(c  -  1)  ki2  .  |-iogi^  ) 


L  W 

wi  kiV 


(1  -  c')  [l  -  c/(l  +  ko2)] 
2(c  -  1)(1  -  c'/c) 


dk(k  +  kj)  e 


^  7-^(Tc  -  Tc')  f  c  /,  k+1  \ 

Jrjgl  +  ks'  c 

lcl2  ,  c'  /  k+1  Tl 


-(log|^-ni)-l 


The  second  solution  differs  in  having  as  an  asymptotic  solution  in  the  tamper  a  growing  exponen¬ 
tial  (growing  for  increasing  negative  x),  e~kix.  The  core  solution  is  again  sinusoidal,  differing  only  in 
phase  from  the  first  solution.  Thus,  the  left  contour  must  still  lie  to  the  right  of  the  roots  of  P(k)  at  ± 
iko.  The  tamper  solution,  g(x),  is  to  grow  as  e-^l^.  Thus  G(k)  must  have  a  pole  at  -k^.  (It  may  also 
have  a  pole  at  +kp  the  corresponding  asymptotic  g(x),  e^l^,  will  be  dominated  by  the  growing  expo¬ 
nential.)  To  give  G(k)  a  pole  at  -k^  the  right  contour  must  pass  to  the  left  of  the  pole  of  P(k)  at  -k^. 
Since  the  left-contour  must  always  be  to  the  left  of  the  right  contour,  the  two  contours  must  be  taken 
as  in  Figure  7.  (Other  contour  arrangements  are  possible,  e.g.,  but  the  solutions  so  obtained  may  be 
represented  as  linear  combinations  of  the  two  solutions  obtained  from  the  contours  of  Figure  6  and  7, 
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Figure  13. 


Deforming  the  contours  of  Figure  7  so  as  to  permit  simplification  of  the  integrals  gives 
this  form: 


Figure  14. 


Taking  as  before: 


log  Pjk)  =  -  ^  /  log  P(k')  =  log  F(k)  +  constant 

Li 

1  r  dk' 

log  PR(k)  =  ^  J  JTTk  constant 


The  integral,  log  PR(k),  may  be  broken  up  into  pieces  which  have  been  evaluated  previously. 

/  V  1  r  dk'  .  /  c  ,  1  +  k'  f  dk'  /  _c^  1  +  k’"^ 

log  PR(k)  "2,11  J  k '  -  k  \  2k'  1  -  k'  ‘  7  ’  2ni  J  k'  -  k  V  2k'  ^  1  -  k'  / 


1  ds 

Jo  s(l  -  ks)  ^0 


1 

2jri 


decaying 


log 


1  4-k^\ 
1  -  kV 
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The  last  term  has  been  evaliiated  in  getting  log  PR(k)  for  the  decaying  solution. 


I  r,  \  _  1  r  ds  „  ds  ds 

og  r(  )  ^  TT  J  s(l  -  ks)  J  s(l  -  ks)  J  s(l  -  ks) 


IT  s(l  -  ks) 

=  ^  Jt^'^c-Tc')  +  2/  ^ -log(k2-ki2). 

•^o  o  o 

log  G(k)  =  i  r^s  ■  "^c' )  ■  -  k2)  +  log  B'  (2,24) 

B  ^ 

It  may  be  observed  that  the  G(k)  here  obtained  differs  by  a  factor  of  ^  ^  from  the  G(k)  previ¬ 
ously  obtained.  Since  the  ratio  of  F(k)  to  G(k)  is  the  same,  the  two  F(k)^s  must  differ  by  the  same 
factor.  We  may  therefore  write  log  F(k)  immediately 


log  F(k) 


k  ds  ,Ti  \  1  B'kp^ (1  -  cQ 

»r  Jo  1  +  ks  ■  ^c')  +  ki2(k2  +  k^'^Xc  -  1) 


B'  is  again  to  be  evaluated  to  give  the  asymptotic  sine  solution  unit  amplitude. 

f(x)  =  sin  kQ(x  -I-  Xj^)  -i-  h(x),  x  >  o,  h(x)--o  as  x— ^  ® 


1  L^l 


21  Ik  -  iko  k  +  ikg 


+  H(k). 


lim  [log  F(iko  +  «  )  -  log  F(-  ik^^  +  « )]  =  log  (-  1)  +  2ikoXj, 


2Lkn  n  ds 


^ikp  r 

^  Jo 


(Tp  -T.O  +  log  (-  1) 


1  A _ ^ 

^  Jn  1 


(Tc  -  TcO  (xi<  0  since  T  <  for  0<  s<  1) 


lim  [log  F  (ik^  -h  €  )  -^  log  F(-  ikQ  -^  € )  +  2  log  €  J  =  -  2  log  2 

,  ki2(c  -  1)  1  ,  _  2(c  -  l)ki2(l  -  c'/c) _ 

=  ^°®ko(l  -  c')  ■  2  [1  -  c/(l  +  ko2)]  (ki2  +  ko2)(l  -  C') 


(cf.  2.19) 
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1  log  ~  ^ 


2ko2(i  -  c')  (1  -  c'/c) 


■  n  I  1  + 

»  r» 


ds  krk^d  -  C ') 

log  F{k)  =  5  I  T-g  (Tc  -  Tc' )  +  log  B'  +  log  ;^-2-(#+L2)(c-1) 


k  /  ds  .  1  ko2(l  -  cQ  [1  -  c/(l  +  ko2)]  (kj^  +  kp^) 

irjl  +  ks'^  ^  2^  2kl2(c  -  l)(k2  +  ko2)2(l  -  c'/c) 


H(k) 


.  .  _  knVki2  +  ko^  .1(1  -  c')[l  -  c/(l  +  kn2)]  rr  1  +  ks  (Tg  *  Tc -) 
■  kj  (k2  +  ko2)  V  2(c  -  1)(1  -  c  '/c) 

kp  Vkl2  +  kn2  /(I  -  c')[l  -  c/(l  +  kn^)]  .  ^  Jp  1 
+  ko2)  V 


ds 


(Tp  -  T.,) 


+  ks^"c  -c 


ki  (k2 


2(c  -  1) (1  -  c'/c) 


k  sin  kpxi  +  ko  cos  ko  xi 
k2Tkn2 


ioo  +  6 


h(x)  =  ^  J  dk  H(k)  ^  j  dk  F(k)  e*“,  {cf.  Fig  12), 
“100+6  -L" 

since  H(k)  is  regular  at  ±  iko  and  F(k)  -  F(k)  -  H(k)  is  single-valued  across  the  — i--  1  cut. 


=  5JrJ 


fdk  ■  ^^^1 

-1 

2 

fpl'-'ks^'^C-Tc-) 

'-si 

log 

1  +  k 

1  -  k 

-  TTi^ 

-00 

l(ki2  - 

k2)(l  -  c') 

1  +  k 

1  -  k 

“ori^  -  1 

“si 

|10g 

|l  +  k| 

ll 

-kl 

1  +  k  1 

\ 

+  TTi  -  1 
/ 

1-k 

where 


D  = 


kpVki^  +  ki^y^d  -  c')  [l  -  c/(l  +  kn2)] 


2(c  -  1)  (1  -  c'/c) 


k  A  ds  ^  . 

■»f/pl  +  ks^^o-Tc') 


h(x) 


_  kicVki^  +  kp2./[l  -  c/(l  +  kn2)](c  -  1)  (1  -  c'/c)  C  kdk  e  _ 

2ko  V  2(1  -  c')  J  k2  -  ki2  7c~k7T7^p^j^ 


-kx 


/c,  k+1 

(,2^°®in 

log  G{k)  =  I  r  (Tp  -  Tp'  )  -  log  (kj2  -  k2)  +  log  B' 

J  r\ 
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rrw  -  ltlVki2  +  ko2  Ac  -  1)  [l  -  c/(l  +  ko^)]~  ^ 1 

'  '  "  kQ(ki2  -  k^)  y  2(1  -  c')  (1  -  c'/c) 


>  ki2  -  k2 


ffjt  /-\s 


G(k)  has  simple  poles  at  ±  ki  and  a  branch  point  at  -1.  We  will  therefore  be  able  to  write  g(x)  as 
g(x)  =  Ae'^1^  +  Be^^^  +  j(x),  j(x)  =  0(ex)  as 


ly,  _ s 

TT  Jo  1  + 


(Tc  “  TcO 


T,  C  tr 

B  =  +  —  e 

2ki 


ij(Tc-Tc.) 


_  k  ds  „  V 

jLi _  «=’' 


c[l  "  c/(2  +  kp^)]  kiV  ki^  +  ko^ 

ki^  -  k^  ko 


2(1  -  c')  (c  -  c') 


k  sinh  kix2  +  kj  cosh  kl 
y  c'  [c  V(1  -  ki2)-2] 


.  .  C  tr 
R(x)=^e 


^  1  -  ^252  -  Tc ')  sinh  ki(x  +  X2)  +  j(x), 


where 


=  i  /  1  -  ki2s2  (’’c  -  ').  (*1<  *2  <  0) 


ki2  s  ds  ,  .  1  iQ.-  (kl^  +  kn2)(c  -  1)  [c  -/  (1  -  ki2)  - 

ffj  rn^^’^c-Tc')  --2^®S  ko2(c/c'- 1)2(1 -cO 


(cf.  2.21) 


eW  =]/%'■ ^  *2)  +  Kx) 


k 

ir  I  1  - 

J  n 


.  1  r  ..  kx  c  ff 

3(X)  =  —  I  dke  e  Jo 

•'-i« 


5(Tc-Te-)__^__3_ 

-k-kj  -k  +  kj 
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-  kx  -  (T  T  ^  ko2(ki2  -  k2)(l  -  c')  2k _ l_j^k _ 

C  f  dke'^  wJo  2  +  ks  (k2  +  ko2)(c  -  l)ki2  ,  c '  i+k 

Zni  ki2  -  k2  ®  ^  2k  ^  1  -  k 


r“  kx  f 

(c  -  c^)  C  kn^(l  ■  cQ  I  kdk  e  e  ^ _ 


kncVki^  +  kp^  /u  -  cQ  (1  "  c  Vc)  ,1  -  c/(l  +  kn^). 

2ki  V  2(c  -  1) 


k  dk  e  o 


e^,  (x<  0) 


We  now  have  two  solutions  whose  asymptotic  forms  are: 


Sin  k„(x  .  XI  .  f  tan-1 


ki(x  +  xg) 


*^1  v^j2  +  ko2  Vc'[c'/(1  -  -  l] 


(cf.  equations  2.17,  2.18,  2.22) 


sin  ko(x  ^  -  kX-\]  «1^Mx  +  X2) 


(cf.  equations  2.25,  2.26,  2.27) 


We  introduce  the  notation, 


j9  =  Vc[l  -  c/(l  +  ko^)] 


^  v/c'[c'/(l  -ki2)  -  l] 


“V^  sm  ko^x  -  XI  .  -  tan  1  ^ 


sin  ko{x  +  XI )  sinh  ki(x  -i-  X2] 


ki2  +  kn^ 


n  (x)  is  ^  ^  Q-  times  the  ‘^decaying  solution’^  first  obtained  (2.14  to  2.23).  ni(x)  is  1/fi  times  the 
°  ki  /5 

^‘growing  solution^’  next  obtained  (2.24  to  2.27).  Subtracting  kin|(x)  from  kinQ(x)  gives 

n2(x)  =  kinQ(x)  -  kini(x) 
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[sin  kQ(x  +  x^)  +  cos  ko(x  +  x^) 


yki2  +  ko^ 


-  ^  sin  ko(x  +  xi) 


k  ki 

^  cos  k^Cx  +  cosh  ki(x  +  X2) 


If  we  now  subtract  ni(x)  from  we  get 


-  ni(x)— I  [cos  k^Cx  +  x^)  ‘  ^  -  sin  kQ(x  +  x^)] 


ki2  +  kn^  .  ,  f  1  j.  -1  ko 

-b- — — ii-  sm  kQ  X  +  XI  -  ^  ’ 


^le‘^l(x  +  X2) 

P' 

We  now  have  two  simple  pairs  of  linearly  independent  solutions,  n(x)  and  n2(x);  n(j(x)  and  n3(x).  For  any 
one  of  these  four  solutions,  hence  also  for  any  other  solution  made  from  them  as  linear  combinations, 
the  asymptotic  solutions  on  the  two  sides  and  the  derivatives  of  the  asymptotic  solutions  have  a  con¬ 
stant  ratio  when  evaluated  at  x  =  -xi  and  x  =  -X2  for  the  core  and  tamper  solutions  respectively. 

derivative  of  asymptotic 

asymptotic  core  solution  (x  =  -Xj)  ^  -k^>9*  ^  core  solution  (x  =  -x^) 
asymptotic  tamper  solution  (x  =  -X2)  kj/9  derivative  of  asymptotic 

tamper  solution  (x  =  -X2) 

the  points,  -x^  and  -X2,  are  both  on  the  core  side  of  the  interface,  -X2  being  the  farther  from  the  inter¬ 
face.  This  property  leads  to  the  following  recipe; 

In  each  medium  the  asymptotic  solution  is  one  of  the  family  of  solutions  of  the  equation:  ( i  +  k^)  n(x)  =  0, 

^  v  -  =  c  (k  may  be  either  real  or  imaginary).  Each  of  the  two  asymptotic  solutions  to  b.e  joined  at  an 
tan“lk 

interface  is  examined  at  its  ''fiducial  point'',  distant  A  x  from  the  interface  on  the  side  of  greater  c. 


ds 

1  -  k2s2 

J  0 

Tc  -  Tc' 

(The  Ax  for  each  solution  uses  its  own  k  which  may  be  either  real  or  imaginary.)  The  two  asymptotic 
solutions,  each  at  its  own  fiducial  point,  have  equal  logarithmic  derivatives.  The  magnitudes  of  the  two 
solutions,  evaluated  at  their  fiducial  points,  have  the  same  ratio  as  their  values  of  the  quantity, 


k2 

c  [1  -  c/(l  +  k2)] 


c  [c/{l  -K  2)  -  1] 


(for  K=  ik) 


*  See  Table  3,  which  gives  c  .  A  X. 
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This  recipe  paraphrases  the  connection-formulae  given  above  identifying  the  two  asymptotic  so¬ 
lutions  on  the  two-sides  of  an  interface.  It  differs  from  a  simple  diffusion  theoretic  boundary  con¬ 
dition  connecting  the  asymptotic  solutions  only  in  so  far  as 

1)  Ax  differs  from  0 

2)  ^  differs  from  a  constant 

This  recipe  connects  only  the  asymptotic  solutions.  Detailed  featxires  of  the  solutions  may  be 
gotten  from  Table  1. 

Symbols  used  in  Table  1. 

In  wtamped  solution 
1 

=  1 

O 


1  +  ko2s- 


-ifl 


^  *  tanh“lk 


-=c,/3'=  ^c[c/(l  -ki2).  i],c  <  1. 


In  tamped  (two-medium)  solutions  the  formulae  have  been  written  for  the  case  c  >1,  c'  <  1.  Other  cases 
follow  by  analytic  extensions. 


_ >  « 

tan’lko 


for  (X2<xi<  0) 
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Each  of  the  four  solutions  is  presented  as  an  asymptotic  solution  in  each  medium  (sinusoidal  or  hyper¬ 
bolic)  to  which  is  added  a  discrepancy  term  (h(x)  for  x>0,  j(x)  for  x<0j.  This  discrepancy  term  may 
be  of  either  sign. 

APPENDIX  I 

ACCURACY  OF  TWO-BOUNDARY  APPROXIMATION 

To  estimate  the  error  introduced  by  neglecting  the  interaction  of  two  boundaries  we  determine  the 
effect  of  this  neglect  in  the  imtamped  sphere  problem  as  a  first  order  perturbation.  The  fundamental 
eigenvalue,  c,  of  the  equation, 

n(x)  =  c  dx'  n(x')yE(|x  -  x  | ),  n(-  x)  =  -n(x).  (i) 

J-a 

we  write  as  c  =  Cq/(1  +  e)  +  0(62)^  where  a  = 

The  integral  operator 

J  jlx'  I  E(|x-x'|) 

we  denote  by  A  . 

Write  R  .  R(x)  =  0  for  x  <  -a 
=  1  for  X  >  -a 
L  ^  L(x)  =  0  for  X  >  a 
=  1  for  X  <  a 
Equation  (i)  becomes 

(1  +  €  -  A  RL)  n(x)  =  0,  valid  for  -a  5  x  ^  a 

n(x)  =  no(x)  +  ni(x)  (ii) 

=  *^rW  +  ^o^ 

where  nj^(x)  and  nj^Cx)  are  the  exact  one-boundary  solutions  satisfying 

(1  -  A  R)nj^  =  (1  -  A  L)nL  =  0 
nj^(x)  =  R  sin  k^x  + 

nL(x)  =  L  sin  k^x  + 

(1  +  6  -  A  RL)nx  =  ( A  RL  -  1  -  €)n^  =  ( ARL  -  1)  (n^^  +  n^^^  -  sin  k^x)  -  €n^ 

=  [a  R  -  1  -  A  R(1  -  L)JnR  +[aL  -  1  -  A  L(1  -  R)]  nj^ 

-^A  -  1  +  A(RL  -  l)j  sin  k^^  - 
=  -  A  ^(1  -  L)nj^  +  (1  -  R)nL  +  (RL  -  1)  sin  kQX  j  -  eriQ 


Then 
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=  -  A  ^  (1  -  Dhjj  +  (1  -  R(hL  +  (R  -  RL  +  L  -  RL  +  RL  -  l)sin  k^xj  -  «  n^, 

=  -  A  [  (1  -  L)hR  +  (1  +  R(hL]  -  «no  (iU) 

Since  nj  must  be  finite,  the  right  side  of  (iii)  must  contain  no  component,  n(x),  satisfying  (ii).  Ne¬ 
glecting  terms  of  order  €  2  we  have 

dx  n(x)-|^A  [(1  -  L)h,j  +  (1  -  R)hL]  +  e  Hoj-  =  0 

€  J"  dx  ao^x)  -  ~  J  RL  n(x)  A  J(1  -  L)hjj  +  (1  -  R)hL  j 

=  -  r*  dx  [(1  -  L)h,j  +  (1  -  R)hL]  A  RL  n(x) 

J-oo 

=  -  r  dx  [^(1  -  L)hj^  +  (1  -  R)hLj  n(x)  (iv) 

«/  ->eo 

The  left  term  of  (iv)  is  roughly  2  a.  The  right  term  is  minus  twice  the  integral  of  the  discrepancy 
term,  hj^  (>0)  starting  from  a  point  distant  2a  from  its  boundary,  with  n(x)  beyond  x  =  a.  The  char¬ 
acter  of  n(x)  in  this  region  may  be  determined  by  taking  c  '  =  0  in  the  decaying  two-medium  solution. 
Its  value  at  the  surface  is 


o')  " 


1  -  c/(l  ^  kn2) 
2 


The  right  term  of  (iv)  will  be  approximately  (-2)  x  - — ^ . 
decay-rate,  about  3-4. 

For  a  tamped  sphere  we  -proceed  in  a  similar  way: 


h(2a)  divided  by  their  combined 


|l  +  €  -  A  [rL  +  (1  -  RL)  ^]j-n{x)  =  0 
n  =  no  +  nj  =  nj^  -f  n^  -  sin  koX  +  n^ 

|l-  A  [R+(l-R)^]}nR  =  |l-  A  [L+(l-L)^]jnL=0 
•jl  +  c-  A  RL  +  y|ni=|A[^^  (Hr  +  Hl  -  sin  koX)  -  «  n^ 

= -[  A  [r  +  ( 1  -  R)  ^  ]  -  1}  Hr  +  A  R(1  -  L)  -  l)  Hr 


+  ^  A  L  +  (1  -  L) 


+  (1  -  L)^J-  lj'nL+  A  L(1  -  -  l^nj 

{1  -  A  RL  +  ^]j-sin  k^x  -  €no 


=  -  A  (1  -  L) 


(R  Sin  k^  +  +  &r) 
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-  A  (1  -  R)  V  +  hL  +  gj 

^  ^  (^i^)  ^  *^0*  ■ ‘"o 

=  (1  -  A  )  sin  kox  -  A  |(1  -  L)hR  +  (1  -  R)hL|-  «n^ 

=  -  ^1  -  a|(1  -  L)  hj^  +  (1  -  R)hLj-  €no 

Hence  as  before: 

«  f  “oW  A  |(1  -  L)hR  +  (1  -  R)hLj> 

Estimating  this  integral  in  the  same  way  as  before  gives,  for  example,  for  c  =  2.0,  c'  =  1.0, 


=  (1  -  A  )  sin  ICqX  -  A  |(1  -  L)hj^  +  (1 


Hence  as  before: 


....  1 


.  2  .5  X. 71  X. 003 

* - 2 - 


For  c'  =  1  and  various  values  of  c,  we  obtain  the  estimates: 


c 

6 

%  in  critical  radius 

1.5 

.0002 

.09 

2.0 

.0015 

.53 

2.5 

.003 

1.0 

3.0 

.005 

1.3 

oo 

.02 

2.0 

The  chief  factor  making  these  errors  small  is  the  rapid  decay  of  h(x).  Taking  the  untamped-solution 
values  as  typical  (they  will  actually  be  somewhat  too  large)  it  would  appear  that  e  will  exceed  .01  only 
for  core  diameters  or  tamper  thicknesses  considerably  less  than  one  mean  free  path. 

Comparison  with  variation  theory  results  gives  about  0.3  as  the  limiting  thickness  for  1  per  cent 
accuracy,  (cf.  Comparison  of  variation  theory  and  end  point  results  for  tamped  spheres,  LADC  -  77) 


APPENDIX  II 

SOLUTION  OF  THE  INHOMOGENEOUS  WIENER-HOPF  EQUATION 

The  Wiener>Hopf  technique  was  shown  by  E.  Reissner  (Journal  of  Mathematics  and  Physics,  Vol.XX 
(1941),  pp  219-223)  to  permit  extension  to  the  inhomogeneous  problem.  We  here  treat  only  the  one  me¬ 
dium  problem  with  the  inhomogeneous  term  confined  to  x  ^  0,  The  extension  to  the  two-medium  prob¬ 
lem  with  an  unrestricted  inhomogeneous  term  is  immediately  obvious.  The  equation  we  wish  to  solve  is: 


(  dx'  n(x' )  K(x  -  x’ )  +  f|(x) 
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where  fi(x)  is  known  and  vanishes  for  x  <  0.  The  Laplace  transform  of  (a),  with  the  notation  used 
previously  is, 

G(k)  =  F(k)  (K(k)  -  1)  +  Fi(k)  =  F(k)  P{k)  +  F^ik),  (b) 


Fj(k)  =  /*  dx  fi(x) 


The  solution  of  the  corresponding  homogeneous  equation  will  be  denoted  by  a  subscript  0. 

Go(k)  =  Fo(k)  P(k) 


P{k)  =  Go{k)/Fo(k) 


We  define  F(k)  such  that 


F(k)  =  Fo(k)  F(k) 

This  introduces  no  singularities  in  F(k)  in  the  right  half-plane  since  FQ(k)  had  no  roots  in  the  right 
half -plane.  Then  (b)  becomes, 

F(k)  P(k)  =  F(k)  Fo(k)  =  i:(k)  Go(k)  =  G(k)  -  Fi(k) 

Thus  -Fi(k)  is  the  right -analytic  component  of  F(k)  Go{k),  which  we  may  write  as 


[F(k)  Go(k)]  R  -  ^  ^ 


where  the  contour  L  lies  to  the  left  of  k  and  of  the  singularities  of  Go(k)  (which  are  entirely  in  the  right 
half-plane)  and  to  the  right  of  the  singularities  of  F(k)  (in  the  left  half -plane). 


[F(k)  Go(k)]R=  -  Fi(k) 


Making  use  of  the  fact  that  .  >  as  well  as  GJk)  is  analytic  in  the  left  half -plane  we  can  show  that 

Go(k) 


equation  c  is  satisfied  by 


[Go(k)F(k)]R  =  -  Go(k)  [Fi(k) 


_ll_  f  G  (k-)  f 

(2»ri)2  j  k'  -  k  ’  J.,  k"  -  k'  Go(k"  ) 

Li  Li 


[Go(k)  F(k)j  ^  Go(k "  \  fi^  ^  ^  k"  -  k  (k '  - 


k  k"  -  k' 
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Displacing  the  contour  L'  to  the  left  of  L"  picks  up  a  residue  at  k'  =  k"  .  The  remaining^ k '  integral 
vanishes  as  it  may  be  displaced  indefinitely  to  the  left,  in  which  direction  the  integrand  decays  as 


=  -  Fi(k)  R  =  -  F^(k) 

The  particular  integral  of  equation  a  has  therefore  the  Laplace  transform 


To  this  may  be  added  any  multiple  of  the  homogeneous  solution,  Fo(k). 

To  extend  this  method  of  solution  to  the  two-medium  problem  requires  only  the  replacement  of 
equation  a  by  the  corresponding  two-medium  equation.  This  leaves  the  form  of  equation  b  and  the  rest 
of  the  solution  unchanged.  To  treat  an  inhomogeneous  term  existing  for  both  x  >0  and  x<0  it  suffices 
to  break  up  the  inhomogeneous  term  into  a  right  and  a  left  side  part  and  treat  each  separately  as  above. 

A  particularly  simple  special  case  of  the  untamped  inhomogeneous  equation  is  that  of  the  albedo 
problem — 
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Then 

r Fi(k)]  f  _ 1 

[Go(k)J  '  2iTi  k'  -  k  (k'  +  a)Go(k') 

1  J_  r  _ dk^^ _ 

Go(-oi)(k  +  oe)  2m  J  (k'  -  k)(k'  +oi)Go(k') 

L' 

In  the  second  term  the  contour  L  may  be  displaced  indefinitely  to  the  left.  Its  integrand  may  be 
written  as 

Const.  \ 

—  *“lk'2) 

Thus  the  k-dependent  part  of  the  integral  vanishes.  The  constant  part  represents  an  admixture  of  the 
homogeneous  solution  to  Fi(k)  and  therefore  may  be  disregarded.  The  general  solution  is  therefore 


In  an  albedo  problem  c  will  be  ^  1  and  A  should  be  chosen  to  make  n(x)  finite  for  all  x,  hence  F(k) 
regular  at  k  =  +ki,  despite  the  pole  of  Fo(k).  Thus 

A-- _ I _ 

Go(-a)(ki+a) 


F(k)  = 


(k  +  ce)(ki  +ci)Go(-a 


The  density  of  emergent  neutrons  in  the  albedo  problem  as  a  function  of  ju,  the  cosine  of  the  angle  of 
emergence,  is 

N(//)  =  c  f  dx  n(x)e 
•^o 


and  is  therefore  given  directly  by  the  solution  F(k). 


^  tW-*W 

Solution  Asymptotic 


CORF  VjuOTION 


:in  K,(x+x,) 


Q(x)-^(x)  'n(o)  =  {(o).^(o)  -H(o)^-p(x)dz 

i#;  [  '•'• 
...J-J-l) 


(<|  Si-nf  K,(xi-x,) 


-i  sin 

V 

=  i  Sin  kJxtJt.) 
P 

f^CoSIt.fxM,) 


ic/c-0  ’^'^0 


j  Sin^  K,X, 


,  -|{,(xfxj 


^  c«,  K,ix*l,) 

P 

=  tf.(».-h,) 


(ntjofive .  <  f  H(o)>o) 


/?'  I  fj 

— !^  [JkO^fc 

p  '  O 


_S1/u£!C 

H(o)*.^L'<.l2(C'')(«A  ttI 


I  [  Sw>  Ml 


^Cej^  K|(z+Xi) 


-KaI 

«>  l|a/c-.)(c-t')  trJ. 


-I  r  ^  r  (i-t'] 

K.[  >^1  Hc->)(f 


Cos  KoX. 


AECD  -  2056 


TA 


0 

LENfetH  OF  Bi: 
IN  toFe 


j(o)  = 

-•p 

MCA  OF  OKCFCWCY 
IN  TAHFfF 

J'(*) . 

0>Wx 

FBMCirNOTH  -OF 
PIStFf  PANf’'  w  T»r«« 

!.(*) 

PreCNEPANCY  (nsoaTivE)  IH  6»«E 
(xio) 

jM 

PiStPEPAHCY  (PbSlTNf)  m  TAMFflf 
(ZSO)  j 

,  ,  , 

Si(i-c'Xc-c')  (C-i)  “ 


(c  - 1  i 


[t-iXc-c')  k4 e 


I  Sl’-'W 


-  K,3t»  j 


fC'i; 


2 


Table 


Table  2,  (continued) 


Table  3.  eilxic.c*) 
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